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Abstract 

It is demonstrated that the properties of light stored in a four-level atomic 
system can be modified by an additional control interaction present during 
the storage stage. By choosing the pulse area of this interaction one can 
in particular continuously switch between two channels into which light is 
released. 
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In recent years a number of nonlinear optical phenomena in weak fields have been inten- 
sively investigated. They are connected with a coherent excitation of an atomic medium, the 
optical properties of which can then be drastically modified. An important example is the 
electromagnetically induced transparency []]], possibly additionally controlled in time, which 
is manifested as light slow-down or even its storage and a controlled release . The 

simplest realization of those effects occurs in atomic systems with three active states in a A 
configuration. Extending such a configuration by adding a coherently coupled fourth level 
opens new possibilities of an external control of such processes [3,p|, |10| , |lll . In our previous 
paper we have shown that in a double A system it is possible to change the light fre- 
quency of the stored light or even to release two pulses of different frequencies by applying 
two control fields, properly chosen and delayed in time. 

Light storing in the form of atomic coherences joins the advantages of the efficiency of 
light as an information carrier and of an atomic medium as an information store. Thus 
it might be a question of a practical importance how to modify in a controlled way the 
properties of the released light by processing the atomic medium during the storage stage. In 
the present work we investigate two possibilities of controlling the released pulse or pulses by 
modifying the atomic coherence due to the stopped light. In the case (a) the lower, initially 
empty state of a typical A system is additionally coupled to a fourth state by another laser. 
In the case (b) the initially occupied state of the A system is coupled by some kind of an 
effective interaction (e.g., magnetic or two-photon electric coupling) with another state of 
the same parity. We show that the Rabi oscillations due to the new interactions modify in 
a coherent way the properties of the released light. 

We consider a quasi one- dimensional medium of four-level atoms with three lower 
metastable states \b > and |c > and \d > and an upper state \a > (Fig.[I]). The posi- 
tion of an atom is described by the variable z, which is considered continuous. The states 
b, a and c constitute a typical A system with the weak signal field 1 and a strong con- 
trol field 2. In the case (a) the state c is additionally coupled with another state d by a 
laser field 4. In the case (b) the state a is coupled with a fourth state d by a weak signal 
field 3 while the states b and d are coupled by some effective coupling 4. The interaction 
Hamiltonian in the case (a) is V = —dJ2j=i.2,A e j cos 4 ) j^ with <pj = uijt — kjz, €j = €j(z,t) 
being slowly varying envelopes and d - the dipole moment operator; we have assumed that 
all the fields have the same linear polarization. In the case (b) the Hamiltonian reads 
V = —dJ2j=i,2,3 e j cos 4 ) j + l(iU)\b >< d\ — {iU)\d >< 6|]cos0 4 , where for numerical rea- 
sons we have made a physically insignificant assumption that the matrix element ill of the 
effective interaction 4 is imaginary. The matrix elements of the dipole moment d\ = (d) a b, 
d2 = (d) ac , <i 3 = (d) a d, c?4 = (d) cd are taken real. Resonant conditions concerning all the 
couplings are assumed. 

The evolution equation ihp = [H, p] for the density matrix p = p(z, t) for an atom 
at position z, after making the rotating-wave approximation, transforming-off the rapidly 
oscillating factors: p ab = a ab exp(-20i), p ac = a ac exp(-i0 2 ), Pbc = cr bc exp[z(0 a - 2 )], 
Pdb = cr d6 exp(-i0 4 ), p dc = a dc exp[i((p2 - fo)}, Pad = o ad exp[z(</> 3 )], p i{ = an, and after 
adding relaxation terms describing the spontaneous emission within the system, takes the 
form in the more complicated case (b) 

io aa = --rr£idi(<7 ba - a ab ) + -Vt 2 ((y ca - a ac ) - —e 3 d 3 (a da - a ad ) - i(T^ + Y a c + T a d )o aa , 
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i&bb = -^ldl((T ab - (T ba ) + ^ iU ( a bd + &db) + i^aa, 
1 

i<j CC = 2^2(o-ac - (yea) + i^c^aa, 

i(?dd = --^e 3 d 3 (a ad — ada) — 2% iU ( abd + a db) + i^d^aa, 

i°ab = -^eidi(a bb - a aa ) + ^ 2 a cb - ^e 3 d 3 a db + ^Ua ad - l -(Tl + T a c + T a d )a ab , (1) 
11 1 % 

i<j a d = --^t\d\O bd + -VL 2 <J c d - -^t?,dz(a dd - a aa ) - —iUa ab - -(F% + r" + r^)(T arf , 

..1 1 1 

^crbc = —-^\d\o ac - -Vt 2 (? ba + 7^iUo- dc , 

itibd = -^t\d\o ad + —e 3 d 3 a ba + t^£7 (ct^ - a bb ), 

i<?cd = 2^20"ad + -^zd z a ca - -iUa cb , 

where T£ is the decay rate of the state \a > to \b >, etc., and ^2 = —e2d 2 /h is the Rabi 
frequency corresponding to the driving field 2. 

The corresponding equations for the case (a) are obtained from the above set (1) by 
interchanging the indices b <-> c, 1 <-> 2, by setting €3 = 0, rf 3 = 0, = 0, by replacing the 
effective coupling iU by = —64^4 (when multiplied by a^c, <r da , cr ac , &db, &bc, &cc and a dd ) and 
by changing the sign of 004.. 

The propagation equations for the signal field 1 (in the case (a)) and for both signal 
fields (1) and (3) (in the case (b)), in the slowly varying envelope approximation and in the 
conditions of the resonance read 



dei lde l u x 

H = -™i — 

oz c at e c 



-K- + —57 = -iMfe ^da, (2) 

c ot e c 

where N is the atom density and eo is the vacuum electric permittivity. Similarly as in earlier 
papers, we have neglected propagation effects for the driving fields, i.e. we take €2,4 = ^2,4(t)- 
Eqs (1) and (2) have been solved numerically in the moving window frame of reference: 
t' — t — z/c,z' — z. Switching the driving field 2 on and/or off was modeled by a hyperbolic 
tangent, while the additional pulse 4 was taken rectangular. The initial probe pulse was 
taken as the sine square shape 

ei(0, t) = ei sin 2 [7r(t — ti)/(t 2 — Ti)]Q(t — ti)©(t 2 — t), while the initial condition for the 
atomic part was a bb (z, 0) = 1, with other matrix elements equal to zero. 

We have performed model computations for data being of realistic orders of magnitude, 
however without making attempt to imitate any real atom. The atomic energies were E a = 
-0.10 a.u., E b =-0.20 a.u. E c =-0.18 a.u. with E d = -0.22 a.u (in the case (a)) and E d = 
E b + 10~ 7 a.u. in the case (b) (the latter value is of order of a magnetic energy splitting). 
The relaxation rates for the spontaneous emission from the level E a to E b , E c and in the case 
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(b) also to Ed were taken equal to 2.4 x 10~ 9 a.u., from which the dipole moments have been 
calculated. The dipole moment for the electric transition c <-> d was taken —2.74 x 10 _1 a.u., 
which corresponded to a negligible width of the level E c . The length of the atomic sample 
was 2.5 x 10 7 a.u. (1.3 mm) in the case (a) and 3 x 10 7 a.u. (1.6 mm) in the case (b) and its 
density 3 x 1(T 13 a.u. (2 x 10 12 cm 3 ). The initial signal pulse length was 10 11 a.u. (2.4 /xs) 
and eio = 10 -10 a.u. (which corresponded to the power density of 3.5xl0 -4 Wcm" 2 ); the 
maximum value of the amplitude of the control field 2 was 1.2 x 10~ 9 a.u. (50 mWcm -2 ). 
The values of the b — d effective coupling U in the case (b) were of order of 10~ 10 a.u. while 
in the case (a) we took a coupling with e 4 = 2 x 10~ 9 a.u. 

In Fig.0 we show the released part of the pulse 1 as a function of the local time t' for 
different values of the area of the pulse 4 (case(a)). The pulse can be lowered, completely 
damped or its sign reversed depending on the final phase of the Rabi oscillations between 
the levels c and d. Of course the final results do not depend on particular time instants 
of switching the interaction 4 on and off, provided that the pulse arrived after the signal 
pulse 1 has been stored and before the release stage has started. The presence of the Rabi 
oscillations becomes clearly visible in the situation in which the control pulses 2 and 4 
partially overlap. In this case, with the latter pulse being now by an order of magnitude 
stronger than before, the restored pulse is constructed of parts freed in those intervals of the 
Rabi period in which the coherence cr foc differs significantly from zero. The Rabi oscillations 
between the levels c and d are thus imposed on the leaving signal pulse (see Fig .[3D. 

The Rabi oscillations due to the additional control field may be used not only to destroy 
in a reversible way the atomic coherence a a b necessary to release the pulse 1 (case (a)) but 
also to create a new coherence a a a which can be converted into a new pulse 3 (case (b)). In 
Fig.^ we show the shapes of the two signal pulses 1 and 3 for different values of the area of 
the control pulse 4, switched on and off in the storage stage. If the multiple of tc 

only the pulse 1 is released, with its sign being changed in the case of an odd multiple. If 
the area is an odd multiple of | only the pulse 3 appears, alternatively with a changed sign. 
For pulse areas being not a particular multiple of | both pulses 1 and 3 are released, their 
heights being under control. 

As in the previous papers the problem can be analyzed in terms of dark state polaritons. 
Such an analysis allows one to describe the whole process of light storing in a single A 
system in terms of a shape preserving solution of the Maxwell-Bloch equations (Eqs (1,2)), 
the components of which, i.e. the signal field and the atomic coherence, adiabatically turn 
one into another. In the case of a four-level system the evolution could not in general be 
fully adiabatic, which means that bright-state polaritons must appear at some stage of the 
process and are later damped [JEJ. Thus the dark-state polaritons at the initial and final 
stages are not identical. 

The approach of Ref. [12[] generalized in our case (b) leads to the following results. 
One can attempt to solve Eqs (1) and (2) perturbatively (as concerns signal fields), in an 
adiabatic and relaxationless approximation. The stage of light stopping occurs as in the 
case of a three-level system: the polariton solution 



describes an adiabatic conversion of the pulse 1 into the coherence and the sign correction 



fi 2 ei + 



o"fec r d 2 




(3) 
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guarantees that for large \Q 2 \ we get ip = e\ (we assume that e 2 > 0). After pulse stopping, 
say at time instant t\, the control pulse 4 is switched on and is present up to the time instant 
t 2 . As a consequence the density matrix evolves and at t = t 2 we obtain in the case (b) 
a hb = cos 2 9, add = sin 2 9, a^d = — singcosg, a bc = (Jbc(ti) cos9, a dc = — &bc(ti) sing, where 
9 = U ^ t2 2h tl ' ) is the pulse area. 

At time instant £3 (t 3 > t 2 ) the control field 2 is switched on in order to release the 
trapped pulse. In the assumed approximations a^, Odd and a^d do not change any more. 
The pulses 1 and 3 satisfy the equations 

/9 d . 13 1^,, , jS 

{ m +c d-J e ' = a 2 mn 2 ^ M ^ (4) 

where j,k = 1,3 and M n = -^^cos 2 g, M 13 = 2N ^ ld? ' sin 9 cos 9, M 31 = 

™>*&<L sin 9 cos 9, M 33 = -H^l sin 2 9. ° 

Eqs (4) can be decoupled by a linear transformation. One of the solutions can be shown, 
similarly as in previous papers, to be a shape-preserving solution traveling with a time- 
dependent velocity 

V (t) = C 2N{d\u Jl cos 2 e+djuj 3 sin 2 gj " 

(The other solution is zero due to the initial conditions at t — t 2 .) The polariton, being a 
combination of two fields and two coherences, has the form 



* = Jdlui cos 2 9 + diu 3 sin 2 9 



toi) d x 



1 + Jg| {djuj! cos 2 9 + d 2 u 3 sin 2 9)\ d ±\ 



diex cos 9 - d 3 e 3 sm9 2N . 

x 32 9/1 , M ^3 H FT i^bc cos 9 - a dc sm9) . (6) 

rffcj! cos 2 g + 6^3 sin 2 g e fi 2 

The solution (6) has been normalized so that it is equal to the solution (3) at t = t 2 , 
i.e. tf(t 2 ) = ^^[(rUh) cos 9 — Gdc (£2) sing]. However, the final form of the polariton 
is, again for large e 2 , a combination of the signal fields with the coefficients different from 
cosg and sing, except in the case of di = <i 3 and ui = u 3 . In particular for 9 = — | one 

finds that \1/ — > ^^£3 \ d l fd 3 \ instead of 63 ^d 3 \ ■ ^ s described in detail in our previous paper 
||12|| , this means that the evolution cannot be fully adiabatic and bright state polaritons 
(which are later damped) must be invoked. In the case (a) the nonadiabatic element of the 
evolution is even more conspicuous: only the part of the transformed coherence, namely that 
proportional to cosg is active in the release stage and turns adiabatically, being a shape- 
preserving solution, into cosgei. The other part of the excitation, namely that proportional 
to sing, "survives" the release stage inside the medium unless relaxations in the states b and 
d are taken into account. 

In summary, we have demonstrated that a modification of the atomic coherence due to a 
stopped light pulse can be used as a new way of changing the properties of the released light. 
One can in particular release two pulses of different frequencies or polarizations, with their 
envelopes being regulated in a continuous way. This may serve as a kind of a switch which 
allows one e.g., to continuously steer the information by sending it to particular channels or 
to temporarily hide it. 
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FIGURES 



FIG. 1. Level and coupling schemes; the indices 1 and 3 refer to signal fields and 2 and 4 - to 
control fields. 

FIG. 2. The field amplitude of the released pulse 1 as a function of the local time t' for different 
pulse areas of the control field 4 in the case (a): curve 1: 0, curve 2: |, curve 3: |, curve 4: |, 
curve 5: §, curve 6: curve 7: tt, 

FIG. 3. The field amplitude of the released pulse 1 as a function of the local time t' in the case 
(a) for overlapping pulses 2 and 4: curve 1: part of the pulse transmitted before light storing, curve 
la: the released pulse in the absence of the additional coupling field, curve lb: the released pulse 
in the presence of the additional coupling field, curve 2: the control field 2, curve 4: the additional 
coupling field 4. The values of the fields 2 and 4 have been reduced by the factors of 40 and 20, 
respectively. 

FIG. 4. The amplitudes (in 10~ n a.u.) of the signal fields 1: full line, and 3: dashed line, as 
functions of the local time t' (in 10 11 a.u.) in the case (b) for different values of the pulse area of 
the interaction 4. 
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